arXiV: astro-ph/0512610v 1 26 Dec 2005 


Massive Star Birth: A Crossroads of Astrophysics 

Proceedings lAU Symposium No. 227, 2005 (c) 2005 International Astronomical Union 

R. Cesaroni, E. Churchwell, M. Felli & C.M. Walmsley, eds. DOI: OO.OOOO/XOOOOOOOOOOOOOOOX 


Physical and chemical structure of dense 
cores in regions of high mass star formation 

Igor Zinchenko^, Lev Pirogov^, Paola Caselli^, 

Lars E.B. Johansson^, Sergey Malafeev^ 
and Barry Turner^ 

^Institute of Applied Physics, Russian Academy of Sciences, 

Ulyanov str. 46, Nizhny Novgorod 603950, Russia 
email: zin@appl.sci-nnov.ru, pirogov@appl.sci-nnov.ru, malafeev_s@mail.ru 
^Arcetri Astrophysical Observatory, Largo E. Fermi, 5, 50125 Firenze, Italy 
email: caselli@arcetri.astro.it 

^Onsala Space Observatory, S-43992, Onsala, Sweden 
email: leb@oso.chalmers.se 

^National Radio Astronomy Observatory, 520 Edgemont Road, 

Charlottesville, VA 22903-2475, USA 
email: bturner@nrao.edu 

Abstract. We found that in regions of high mass star formation the CS emission correlates well 
with the dust continuum emission and is therefore a good tracer of the total mass while the N 2 H'*’ 
distribution is frequently very different. This is opposite to their typical behavior in low-mass 
cores. The behavior of other high density tracers varies from source to source but most of them 
are closer to CS. Radial density profiles in massive cores are fitted by power laws with indices 
about —1.6, as derived from the dust continuum emission. The radial temperature dependence 
on intermediate scales is close to the theoretically expected one for a centrally heated optically 
thin cloud. The velocity dispersion either remains constant or decreases from the core center 
to the edge. Several cores including those without known embedded IR sources show signs of 
infall motions. They can represent the earliest phases of massive protostars. There are implicit 
arguments in favor of small-scale dumpiness in the cores. 
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1. Introduction 

Many important aspects of high mass star formation (HMSF) remain unknown or 
poorly known. This is caused by the much more complicated appearance of HMSF regions 
in comparison with their low-mass counterparts and by the relative difficulties of their 
studies due to larger distances, etc. Therefore, detailed investigations of the structure 
and kinematics of HMSF cores are of great importance. For example, different models 
of star formation predict different density profiles in the cores and different frequency of 
outflow occurrence. 

An important question in this respect is a selection of reliable tracers of physical 
parameters in HMSF regions. Chemical effects can seriously distort their appearance. 
It is now well established that the central parts of dense low mass cloud cores suffer 
strong depletion of many molecules (in particular, CO, HCO+ and CS) onto dust grains 
(Caselli, Walmsley, Tafalla, et al. 1999; Caselli, Walmsley, Zucconi, et al. 2002; Kramer, 
Alves, Lada, et al. 1999; Willacy, Langer, & Velusamy 1998; Jessop & Ward-Thompson 
2001; Tafalla, Myers, Caselli, et al. (2002); Bergin, Ciardi, Lada, et al. 2001), suggesting 
that CS (so far considered a high density tracer) does not actually probe the central 
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core regions. On the other hand, N 2 H''' is usually an excellent tracer of dust continuum 
emission (Caselli et al. 2002), implying that this species depletes at higher densities 
compared to CO, despite of the same binding energies (Oberg, van Broekhuizen, Fraser, 
et al. 2005). 

The situation is less clear in HMSF cores where the mechanisms of chemical differenti¬ 
ation are probably different (e.g. not linked to molecular freeze-out). There were several 
studies which show chemical variations in a few similar regions (Ungerechts, Bergin, Gold¬ 
smith, et al. (1997); Bergin, Goldsmith, Snell, et al. 1997). Here, we present observations 
of a statistically significant sample of HMSF cores. 


2. Chemical differentiation in HMSF regions 

A comparison of our CS (Zinchenko, Mattila & Toriseva 1995; Zinchenko, Pirogov & 
Toriseva 1998; Pirogov et al., in preparation) and N 2 H+ (Pirogov, Zinchenko, Caselli, 
et al. 2003) maps in many cases shows striking differences between them. Our further 
studies of dust continuum emission and other high density tracers (HCN, HNC, HCO+, 
their isotopes, etc.) reveal certain relationships between them. The details of these ob¬ 
servations are presented elsewhere (Pirogov et al., in preparation, Zinchenko et al., in 
preparation). Here we summarize the main results. 

First, we found a good correlation between CS and dust emission in most cases. At 
the same time N 2 H''' distributions are frequently very different. An example of such 
comparison is shown in Fig. ^ The millimeter wave dust (and CS) emission peaks are 
also usually associated with strong FIR sources (in particular, from the IRAS catalogue). 
On the other hand the N 2 H+/CS ratio usually reaches the highest values towards clumps 
without such association. 




Figure 1. Maps of CS(5-4) and N2H^(l-0) emission (contours) overlaid on the 1.2 mm dust 
continuum map (grey scale) in G285.26-0.05. Star marks the position of the IRAS point source. 
All maps are obtained at SEST. 

This behavior is opposite to that typical for cold low-mass cores. The variations of 
the CS/N 2 H+ intensity ratio reach more than an order of magnitude. It is worth noting 
that the N 2 H+ optical depth is usually small (Pirogov et al. 2003) while the CS optical 
depth can be rather high. This means that the saturation effects can only decrease the 
apparent variations of the CS/N 2 H+ ratio. 

The behavior of other high density tracers varies from source to source but most of 
them are closer to CS. An exception is HNC which is an intermediate case in this sense 
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between CS and N 2 H+. In Fig. we present as an example maps of several high density 
tracers in S255 overlaid on the dust continuum map. It shows two continuum peaks 
almost equal in intensity. The nature of these two components is different. The central 
one is associated with a luminous cluster of IR sources, whereas toward the northern 
one an ultracompact H ii region (G192.58-0.04) was detected. CS and HCN follow the 
dust distribution quite well, while N 2 H+ is very different. HCO'*' map in S255 looks 
similar to N 2 H+ but in some other sources it is closer to CS. The same can be said about 
ammonia. A comparison of the C^®0 and dust continuum data for this source shows 
no sign of CO freeze out which suggests that the dust temperature is above 20 K, the 
sublimation temperature of CO. 



Figure 2. Maps of S255 in various molecular lines (contours) obtained with the OSO-20m radio 
telescope overlaid on the map of 1.2 mm dust continuum emission (grayscale) obtained by us 
with the IRAM-30m radio telescope. 


We studied variations of kinetic temperature in these sources by CH 3 C 2 H observations. 
In most cases no significant temperature difference between the CS and N 2 H+ peaks was 
found (Malafeev, Zinchenko, Pirogov, et al. 2005). Typical temperatures which we derived 
for these peaks are ^ 30 — 40 K. E.g. in S255 we found kinetic temperatures of ~ 40 K for 
both peaks mentioned above. Therefore, the observed chemical differentiations cannot be 
explained by molecular freeze-out as in low mass cores. 

Density differences between the CS and N 2 H+ peaks are also not very significant as 
follows from our density estimates based on CS and CH 3 OH LVC modelling. At the same 
time the gas densities at CS peaks are somewhat higher than at the N 2 H+ ones and in 
general there is a correlation between CS emission intensity and gas density. 

One possible explanation for the observed chemical differentiations was proposed by 
Lintott, Viti, Rawlings, et al. (2005). They suggested that the enhancement of CS/N 2 H+ 
abundance ratio may be related to the high dynamical activity in these regions which 
could enhance the rate of collapse of cores above the free-fall rate. Consequently, high 
gas densities would be achieved before freeze-out had removed the molecules responsible 
for the N 2 H+ loss, while the high densities promote CS formation. Our data show that 
in most cases the line widths at the CS peaks are somewhat larger than at the N 2 H+ 
peaks which is consistent with this model although in S255 the reverse is true. 


3. Internal structnre and kinematics of HMSF cores 

3.1. Radial density profiles 

We studied the density profiles in the cores using the molecular line and dust continuum 
maps. From the N 2 H+ data we derived the ~ r~^ density profile assuming constant 
N 2 H+ excitation and abundance along the radius (Pirogov et al. 2003). However, this 
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assumption is not justified and probably more reliable conclusions can be drawn from 
dust observations. The analysis of the dust continuum maps was similar. We fitted the 2D 
intensity distribution for nearly circular clumps with a convolution of a power law func¬ 
tion with a gaussian beam. The details of this analysis are presented elsewhere (Pirogov 
et at, in preparation). An example of this fitting is shown in Fig.|21 
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Figure 3. An example of the radial dependence of dust emission intensity fitted by a power 
law function convolved with a gaussian beam (left panel) and the radial dependence of density 
derived from LVG analysis of the CS J = 2 — 1 and J = 5 — 4 lines (right panel). 


In total, the data for 17 clumps obtained at SEST were analyzed. The power law index 
is p = —1.0 ± 0.1 for clumps with embedded IRAS point sources and p = —0.6 ± 0.1 
for clumps without such sources. Assuming temperature gradient in the former sources 
(described by a power law index of —0.4, see below) and isothermal latter sources, this 
corresponds to ~ density profiles for both types of clumps. 

Such density profiles are in agreement with several other studies (e.g. Beuther, Schilke, 
Menten, et al. 2002; Mueller, Shirley, Evans II, et al. 2002). They rule out some theo¬ 
retical models (e.g. so-called “logatropic” one) but are consistent with the “standard” 
theory of star formation (Shu 1977). 

It is worth noting that these dependencies are obtained for mean volume densities 
which are derived from column densities. Densities found for regions of molecular emission 
from excitation analysis are usually significantly higher and have different dependence 
on radius (as shown in Fig. 0 ). This can be a consequence of small scale dumpiness in 
these regions. The volume filling factor should be ^ 0.4. 


3.2. Radial temperature profiles 

We studied the temperature distribution in the cores on the basis of high quality CH 3 C 2 H 
J = 13 — 12 maps obtained at IRAM 30-m telescope (the beam width at this frequency is 
~ 12"). An example of temperature map (for S140) is shown in Fig.0 These maps could 
be deconvolved into several distinct clumps and the dependence of temperature on the 
projected distance from the clump center was analyzed. An example of this dependence 
for the main clump in S140 is also shown in Fig.0 

These dependencies were compared with model calculations of temperatures which 
would be derived for optically thin CH 3 C 2 H emission from a spherically symmetric cloud 
with a power law temperature and density gradients (CH 3 C 2 H was assumed to be ther- 
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Figure 4. Distribution of kinetic temperature in S140 as derived from the CH 3 C 2 H J = 13 — 12 
data obtained in 2004 with the IRAM-30m radio telescope (left panel) and the dependence of 
this temperature on the projected distance from the clump center for the strongest clump (right 
panel). Positions of IRS 1, 2 and 3 (Evans II, Mundy, Kutner, et al. 1989) are marked. Apparent 
variations of kinetic temperature at the periphery of the map are unreliable and are probably 
caused by low signal to noise ratio. 


malized at densities n > 10^ cin“^). We found that the radial temperature dependence 
in these clumps can be fitted by a power law with indices between -0.3 and -0.4. 

This dependence is in agreement with the theoretically expected one for a centrally 
heated optically thin cloud. Some other recent studies found somewhat different depen¬ 
dencies. E.g. Fontani, Cesaroni, Caselli, et al. (2002) derived steeper temperature gra¬ 
dients from comparison of temperatures found from different species and at the same 
time argued for rather constant temperature in regions of CH 3 C 2 H emission. The results 
presented here do not support these conclusions. 

3.3. Radial velocity dispersion profiles 

Our CS and N 2 H+ data show that the velocity dispersion in HMSF cores either remains 
constant or decreases outwards (Zinchenko 1995; Lapinov, Schilke, Juvela, et al. 1998; 
Pirogov et al. 2003). A plausible explanation could be found in a higher degree of dy¬ 
namical activity of gas in central regions of HMSF cores, including differential rotation, 
infall motions and turbulence due to winds and outflows from massive stars. It is worth 


noting that Caselli & Myers (1995) derived an opposite trend from comparison of line 
widths and sizes of the emission regions for different species in the same source. The 
reason for this discrepancy is not clear yet. 


3.4. Core rotation and associated outflows 

Many HMSF cores mapped in N 2 H+ show systematic velocity gradient fields, implying 
nearly uniform rotation (Pirogov et al. 2003). There is a correlation (cc = 0.9) between 
direction angle of total velocity gradient and elongation angle. This fact justifies an 
assumption that elongation of cores and clumps could be due to rotation. However, the 
ratio of rotational to gravitational energy is only ~ 0.01 on the average, so that rotation 
should not play a significant role in core dynamics. 

The frequency of occurrence of high velocity outflows in the studied sample was found 
to exceed ^ 40% using SO line wings as an indicator (Zinchenko 2002) which is in 
agreement with other available estimates. 
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3.5. Infall motions and massive protostars 

In several cases molecular line profiles (especially HCO'*') show features indicative of infall 
motions (red-shifted self-absorption). This is a case both for clumps with and without 
strong embedded IR sources. However, the latter ones are probably at an earlier stage of 
evolution, and they are characterized by a relatively high N 2 H+/CS ratio, as mentioned 
above. Thus, one can speculate that this ratio can be one of the indicators of the earliest 
phases of massive protostars. 


4. Conclusions 

We have shown that effects of chemical differentiation in warm massive cores are very 
different from those in cold low mass cores. In particular, CS is a good tracer of total 
mass here, while N 2 H+ is not. A relatively high N 2 H+/CS ratio can indicate the earliest 
phases of massive protostars. The internal structure of HMSF clumps is more or less 
consistent with the standard theory of star formation. 
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